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The formalism of Causal Dynamical Triangulations (CDT) attempts to pro- 
vide a non-perturbative regularization of quantum gravity, viewed as an ordinary 
quantum field theory. In two dimensions one can solve the lattice theory analyt- 
ically and the continuum limit is universal, not depending on the details of the 
lattice regularization. 
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1 Introduction 



Two-dimensional quantum gravity has been a fruitful laboratory for studying 
aspects of string theory as well as quantum gravity. One somewhat surprising 
aspect of Euclidean two-dimension quantum gravity coupled to matter in the 
form of a conformal field theory, is that the regularized lattice theory, using the 
so-called dynamical triangulations (DT), can be solved analytically. The details 
of the DT regularization are unimportant for the continuum limit. In fact it 
has been a wonderful example of universality in the Wilsonian sense, the critical 
surface where the continuum limit can be taken being of finite co-dimension in 
an infinite dimensional coupling constant space (see e.g. [1] for a review). The 
lattice regularization known as causal dynamical triangulations (CDT) uses a 
subset of the triangulations used in DT [2, 3]. The original idea was to consider 
a path integral where spacetime histories before rotating to Euclidean signature 
were locally causal, i.e. had non-degenerate light cones (see [4] for a review of the 
CDT approach also in higher dimensions than two). In two dimensions, which 
is the only case we will consider here, the precise relation between the CDT 
triangulations and the DT triangulations was described in [5]. 

There is good evidence of universality of the CDT scaling limit, although one 
does not have the same comprehensive evidence as for the DT case. First, a 
related model, in a certain way more general, the so-called string-bit model [6], 
led to the same scaling limit. Further it was shown in [7] that one could add 
dimers on the "spatial" CDT links without changing the universality class. Thus 
it was somewhat surprising that adding further "dressing", but only along the 
spatial links, seemingly led to new continuum models, depending on a continuous 
parameter ft (to be defined below) [8] . The purpose of this letter is to show that 
also for this general set of models one obtains indeed the standard CDT scaling 
limit. 

2 Defining the model 

The modified CDT model (not to be mistaken for what has later been called 
"generalized CDT" [9]) is most easily defined using a lattice dual to the triangu- 
lation, i.e. a 3 graph with a "time" foliation . Fig. 1 shows the dual CDT lattice 
and its generalization. In this dual picture each vertex represents a triangle in 
the "original" triangulation and each polygon represents a vertex, the order of 
which is equal the number of sides in the polygon. 

In the modified model one allows a dressing of the horizontal links between 
two vertical links by rainbow diagrams. 

Three coupling constants are assigned to the model: to each vertex one asso- 
ciates a coupling constant g, to a vertex with an incident vertical link an additional 
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Figure 1: Modified CDT configuration, dual graph. 



coupling constant h, and finally to each vertex with an incident rainbow link a 
coupling constant 6. The parameter 



governs the density of rainbow links compared to the number of vertical links, i.e. 
"time-like" links in the original CDT-like 3 -graph. In this article we will only 
consider < (3 < 1, which is the range leading to CDT-like theories [8]. 

As shown in [8] one can define and calculate a transfer matrix for this model. 
The result is 



where the index j refers to the number of incoming half-lines which is incident 
from below on the horizontal line at time t and index k refers to the number 
of half-lines leaving the horizontal line at time t. Index k plays the same role 

(2) 

as index j, only at time-slice t + 1. In this way & ik connects outgoing vertical 
half-lines at t to incoming half-lines at t + 1 and Qij incoming half-lines at t to 
incoming half-lines at t + 1. 

0^) is the CDT transfer matrix, already discussed in [2] and analyzed in detail 
in [7]. If 9 = and h = 1 there are no rainbow lines and O*- 2 -* becomes the identity 
matrix and also the CDT transfer matrix. 




(2) 



k 



It is convenient to work with the discrete Laplace transforms of 0, 0^ and 
@( 2 ). To simplify the expressions somewhat we make the following redefinitions 
compared to [8]: 

eg> (2 5 )-^eg ) , Qif (2g)^Q\f. (3) 

The explicit expressions are then: 

e«(x,y) = ^xVeg ) = — T i- T - (4) 



{ ,y) ' (l-x 2 C{x 2 )){l-y 2 C{y 2 )){l- (3~ 2 xy C{x 2 )C{y 2 )) K ) 

&(x,y) = jf^e<V" _1 )e (2 W), (6) 

where the contour encloses cuts and poles and where 

x = 2g6x, y = 2g6y, C(z)= 1 ~^~ —. (7) 
Integrating over the simple pole of one obtains 

Q (X y) = -1 W) 1 (8) 

{ ,y) l-lx{l-x 2 C(x 2 ))(l-y 2 C(y 2 ))l- (3- 2 xyC(x 2 )C{y 2 Y K ) 

where 



* = ( 9 ) 
2-x 



The partition function with open horizontal boundaries after t time steps is 1 
Z(l,k;t) = f(0 (1) (© (2) © {1) )' N ) , (10) 

V / kl 

and the (discrete) Laplace transformed function is denoted Z(x,y) 

Z(x,y;t) = J2^y k Z(l,k;t). (11) 



Lk 



The partition function after t time steps with periodic boundary conditions in 
the time direction is 

Z(t)=tr(e*). (12) 
1 The same continuum limit is obtained by setting Z(l,k;t) = (6* 



'kl 
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3 The continuum limit using the transfer matrix 



As shown in [8] the partition function Z(t) has a singularity at 

Zc = 2ge(0+^) =1. (13) 

We want to take to continuum limit by approaching this singularity. This is done 
in the following way [8]: 

f^+i) =1-^(1^)'. (14) 

The interpretation is that a is the lattice spacing, i.e. the link length in the 
triangulation, and A the cosmological constant, such that the average number 
of triangles is proportional to l/(Aa 2 ). Thus the average "continuum" area is 
proportional to 1/A. 

Until now t has denoted the integer number of time steps in the triangulation. 
We are interested in a limit where we have a finite continuum time T scaling as 

T = ta, (15) 

where a is the lattice spacing defined by (14). We can then write 

Z(T) = tr 0* = tr e~ TH , 6 = e- aH . (16) 

Thus an expansion of 9 to lowest order in a should allow us to determine H. 

If the continuum area is proportional to 1/A we expect the continuum length of 
a time slice to be proportional to 1/(AT). Thus we expect a scaling L oc I a where 
/ is the number of space-like links. We can also enforce this on the boundaries: 

Z(l,k;t)^Z(L ,L T ;T). (17) 

The discrete Laplace transform of Z(x, y; t) has poles in x, y and it is at these poles 
one extracts the continuum function Z(L , Lt;T). These poles are at x c = y c = 1 
for a — > 0. The terms x l and y k in (11) can then be given an interpretation as 
the part of the action coming from a continuum boundary cosmological term 
proportional to X if we scale: 

and thus 

x l ->■ e~ LX for a ->■ 0. (19) 



5 



With this scaling we obtain a relation similar to (17), going from the discretized 
expression to the continuum expression: 

Z(x,y,t) ^ Z(X,Y;T), (20) 

where the continuum analogue of (11) reads 

POD 

Z(X,Y;T)= / dL dL T e- LoX - LTY Z(L ,L T ;T). (21) 
Jo 

We will return to (17) and (20) in the next section. 

We now extract H from G = e~ aH . It is convenient to use the Laplace 
transform (6) of ©. Expanding in a we obtain [8]: 

((1 - aH + 0( ' W ) (st ) = f£- 8 (,, i) 4-H. (22) 

Here ^(u) is the discrete Laplace transform of a function 

^(u) = J2"W)- (23) 
i 

The function 9(x, 1/cu) has a pole in u; at 1 for a — > and it has a branch cut 
located at w 6 [—a;*, a;*], where 

u* = 2^p+^J + 0(a) < 1 for a sufficiently small. (24) 

We can deform the contour to be a small circle around one and an integration 
along the branch cut. The integration around u = 1 allows us to use the expansion 
(18) for x and uj, and we obtain 



dZ 
2tH 



1 | a f\ | ~ t 1 + 3f3 2 )XZ + f3 2 Z 2 



Z -X (Z -X) 2 \ 1 + 



2 



^(Z)+0(a 2 ), 
(25) 

Performing the integration (and ignoring the contribution from the cut) we can 
identify H as 

H(X) = (X 2 -A)^+X, (26) 
and by an inverse Laplace transformation 

H ^ = - L &~M + AL - (27) 

This is precisely the ordinary CDT Hamiltonian, the only difference is that in 
order to obtain it in this form we had to perform a dressing (or renormalization) 
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of the continuum boundary cosmological constant from a value X, corresponding 
to /3 — to the f3 dependent value given in (18). This renormalization of X and a 
similar renormalization of the coupling cosmological coupling constant A in (14) 
is all that is needed to include the effects of the rainbow diagrams. 
The contribution from the cut can be written as 

diuf(x,u)^(u), (28) 

where f(x,u) is integrable in [— u*, w*] and ip{x) analytic in the neighborhood 
of 1 and finite when a — > 0. We cannot view such a function as the Laplace 
transform of any function ^(y/AL) depending on the continuum length L > 0, 
the reason being that the inverse Laplace transformation from (26) to (27) gives 

/ e XZ ^(l - aX) = 5(L)tP(1) - a8\L)fi(l) + ■■■ + 0(a n ). (29) 

Thus we do not associate any continuum physics with the analytic function ip(x) 
defined by (28) 2 . 



4 The Schwinger representation and the contin- 
uum 

In [8] the modified CDT Hamiltonian was not derived using the transfer matrix as 
described above, but rather a so-called Schwinger representation of Z(x, y; t). We 
now show that this method also leads to (27), i.e. the ordinary CDT Hamiltonian. 

The starting point is the following representation of Z(x,y;t) ([8], formula 
(5.19)): 

* / roc \ * _1 

Z(x, y; t) = H / da s e- a ° e ^ a ° x+a ^ J] M90a r ,g9a r+1 ) (30) 

s =0 \ Jo ' r=0 

where 

( t )p {x,y) = Y,h{2x)I k {2y)/p 2k . (31) 

fc>0 

2 Of course a function like ip(uj) would also not contribute to continuum physics if inserted 
in (25). The part of a function ip(u>) defined as in (23) which does contribute to continuum 
physics in (25) is the part which has a continuum Laplace transform, i.e. the part where ip(l) 
in (23) has the form i/>(i/£ — £ c - > il>{\/A-L). Since y/£ — £ c oc a^/A it can at most be the tail 
at infinite I which contributes to continuum physics for a given tp(ui) = w y(0- 
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x and y only appears in the exponential function and we can write 

POO POO 

Z(x,y;t)= da da t e-^ l - x)a °-^ l - y)at F(a ,a t ;t), (32) 
Jo Jo 



where 

't-1 „no \ t-i 



F(a ,a t ;t) = I J" f / da s ) ]Je~^^M^r, 9 ea r+l ). (33) 

\»=1 J ° J r=0 

Since 1 — x oc aX and 1 — y oc aF, (32) states that in the limit where a — > and 
thus Z(x, y; t) ->■ Z(X, Y; T), Z(X, Y; T) is the Laplace transform of F(a , a t ; t), 
t = T/a. Thus, in accordance with (21) we have 

F{a ,a t ;t) <x Z{L , L T ;T), (34) 

where 

1 /l - /5 2 \ 2 T 1 fl-fi 2 ^ 2 



L ° = 2 aa ° \lT]^ ) ' Lt = 2 aat \T+JP ' ' at = T - ( 35 ) 
If we change variables from a s to v? s , 



a \1- 

we obtain 



1 /"°° TT TT \/ ( Pr L Pr+l 1 + f3 2 a r +a r+1 

Z(L ,L T ;T) oc / I d<p s I — e 2 (f)p(g9a r , g0a r+1 ). 

(37) 

The right hand side can be interpreted as a (quantum mechanical) path integral, 
i.e. 

JpWt Z(L , L T ; T) oc (<p \e- TH \<p t ) (38) 

for some Hamiltonian H . We will now proceed to determine H. 
Following [8] we use the notation 

e-^M90ao, gOaJ ~ Up(a , a ± ) e~ s ^ ai \ (39) 

According to [8] 

Sp(a , oti) = ^(«o + «i) ~ 2g0 V («o + /3 2 ai)(«o + /3~ 2 «i) (40) 
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and 

U j3 (a ,a 1 ) 



^4^9 ((a Q + f3 2 a 1 )(a + P' 2 ^)) 1 / 4 
x I 1 -\ ' 



16^v/(a + /3 2 «i)(ao + /3- 2 ai) 



+ 



We now expand in a, with 



A(p = (pi - ip 



(41) 



(42) 



counted as being of order yfa as one has to do in a path integral (here we differ 
from [8]): 



Ay? 2 
~2a 



(3 2 Ay? 4 aA 2 ^. o /2x 
^ + — V^ 2 + 0(a 3 / 2 ). 



l + /3 2 ) 2 2ay? 2 



(43) 



We see that we get a standard kinetic term, justifying Ay? oc -^/a. (Note that the 
Ay? 4 term is not present in 
Similarly, we find 



y/(p r ip r + 1 1 + 

a 1 -/3 s 



^(ao,ai) 



72 



1 / a 

1 + 



ira 



:i+/? 2 ) 2 <pI + ( 7' 

(44) 



(We note that the Ay? 2 term is not present in [8].) 

The Hamilton is finally determined by integrating against a trial state: 



((l-atf)V)(<A)) 



+ 



dy?i 



o y/2na 
2 Ay? 4 aA 



1 + 



(3 2 Ay 2 
1 + /3 2 ; 8y? 2 (1 + /3 2 ) 2 V 2 



(l + /3 2 ) 2 2ay? 2 2 
Carrying out the Gaussian integral, we obtain 



1 + Ay?— + 



d Ay? 2 d 2 



dp 2 dp 2 



iPM- (45) 



H = - 



i a 2 



A 2 

+ ^ - 



1 



n2' 



2<9y? 2 ' 2 

This is precisely the CDT Hamiltonian when changing back to the L variable. 



(46) 



5 Critical arches 

In principle a new behavior could be possible for j3 — > 1 from below, since in this 
case the rescaling of lengths and boundary cosmological constants, as defined by 
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eqs. (18), diverges and it is precisely the limit where the cut will merge with the 
pole in the expression (8) for 6. Let us investigate this case by assuming 

= 1- a r 'B, (47) 

where B is a new physical constant with mass dimension rj. To understand the 
analytic structure of © for a — > 0, i.e. (3 — > 1 from below, we expand the argument 
of the square root related to the Catalan number in the expression for 0: 

Vl - 4£ 2 = aPB(l + aX + \a 2 {K - X 2 ) + 0(a r 'B) + 0(a 3 )) (48) 

From this expression it is clear that that the cut has disappeared from the ex- 
pression even though it hits the pole when expressed in terms of unrenormalized 
variables. To find the Hamiltonian we use the same approach as in Sec. 3, eqs. 
(22) and (25) and write 

§{x) = (1 - a"H + ■ ■ ■ Mx) := a ^-l^- Q( x , (49) 



where v is determined by the expansion, We find: 
dZ 



MX) 



2m 



aPB/2 a A + \{X 2 - AXZ + Z 2 ) 



HZ)- (50) 



Z-X (Z-X) 2 

Thus, if i] > 1 we obtain the same results as before (eq. (25) with = 1) and if 
rj < 1 we obtain a trivial Hamiltonian. rj = 1 just adds the positive constant B/2 
to the CDT Hamiltonian (26). So far we have ignored the contributions from the 
cut. However, arguments like the ones used in Sec. 3 show that the cut will not 
contribute in the scaling limit. 



6 Discussion 

We have shown that the CDT scaling limit is quite universal and independent of 
details of the lattice regularization, as long as we maintain a reasonable "memory" 
of the underlying assumed time foliation. Dressing the spatial slices with a few 
outgrowths should not alter the scaling limit and this is indeed what we have 
proven to be the case. Potentially there could have been a different behavior in the 
limit /3 — > 1 where the rainbow diagrams become critical, but explicit calculations 
showed that it was not the case. The CDT model provides us with a regularized of 
a theory of fluctuating spacetime which is invariant under spatial diffeomorphisms 
and which allows for a time foliation. The simplest such continuum model is a 
Horava-Lifshitz gravity model in two-dimensions where we only keep terms with 
at most second order derivatives of the metric, and one can indeed show that such 
a model has a classical CDT Hamiltonian which when quantized is compatible 
with the H(L) considered in this paper [10]. 
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